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Abstract 

We study the decay B —)■ Xg'j in Randall-Sundrum models with an IR-localised 
bulk Higgs. The two models under consideration are a minimal model as well as 
a model with a custodial protection mechanism. We include the effects of tree- 
and one-loop diagrams involving 5D gluon and Higgs exchanges as well as QCD 
corrections arising from the evolution from the Kaluza-Klein scale to the typical 
scale of the decay. We find the RS corrections to the branching fraction can be 
sizeable for large Yukawas and moderate KK scales T ; for small Yukawas the RS 
contribution is small enough to be invisible in current experimental data. 



1 Introduction 


One of the best studied processes in flavour physics is the inclusive radiative B —)■ 
decay. On the experimental side numerous experiment s (l} j^ provide an ever increasing 
amount of data; leading to the current HFAG average |9] of 


Br(5 ^ 


E.^>1.6GeV 


= (342 ± 21 ± 7) X 10“® 


( 1 ) 


where all contributing experimental results were converted as to correspond to a lower 
photon energy cut of 1.6 GeV. A further improvement of this number can be anticipated: 
the Belle II experiment is expected to be able to measure the branching fraction with 
an uncertainty of about 6% 
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On the theory side, the fact that the rare radiative decay provides both powerful 
check for the Standard Model (SM) of particle physics and is sensitive to physics beyond 
the SM (BSM) fuelled a tremendous effort (see e.g. 11-14 and references therein) to 
understand the intricacies of the 6 —>■ sy transition. The most recent result is given 
by 


Br(S ^ X, 7 ) 


th 


= (336 ± 23) X 10 


-6 


( 2 ) 


It is in very good agreement with experiment, cp. ([^, and therefore provides non-trivial 
constraints to any New Physics model that can generate additional flavour-changing 
neutral currents (FGNGs). 

Extra-dimensional models of the Randall-Sundrum (RS) type are known to have 
a particularly rich flavour phenomenology and can, despite an inherent protection mech- 
give rise to sizeable FGNGs. The characteristic Eve-dimensional metric of 
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anism 

RS models can be written as 


ds^ = {r]^ydx^dx'^ — dz'^) , 


( 3 ) 


in conformal coordinates. Here k = 2.44 ■ 10^® GeV is of order of the Planck scale Mpi. 
The hfth coordinate x is restricted to the interval [1/fc, 1/T]- The boundaries z = 1/k 
and z = 1/T are typically referred to as Planck and IR brane respectively. The a priori 
arbitrary scale T is assumed of the order of a TeV in order to alleviate gauge-gravity 
hierarchy issues 
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One of the main reasons for the popularity of these models is the interplay of (SM) 
flavour and properties of 5D wave functions 19 21 . In particular, mass and GKM 


hierarchies can be related to the strength of the Planck or IR brane localisation of the 
corresponding KK zero-mode wave functions 22 . This intimate relationship of geometry 


and flavour makes the study of flavour physics observables all the more intriguing. For 
most processes like meson mixing [^^4] or electroweak pseudo-observables 25 -27 the 


RS contribution arises (to leading order) from tree-level corrections to dimension-six 
operators, e.g., four-quark operators in the case of meson mixing. 
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In the last few years loop-induced processes, that is processes that to leading order do 
not receive contributions from tree-level diagrams in RS models, have been studied quite 
extensively. Observables that have been investigated include /r —>■ ey (28| - [30| , {g — 2)^ 
Higgs production and decay 


31 32 


33 38 as well as c —)■ uy and c ^ ug 1^. The 


latter process just as fi ^ e'y receives contributions from Kaluza-Klein (KK) states of the 
Higgs (in models where these are present). The subtleties involving their determination 
have only recently been pointed out (40) . 

The decay B —)■ X^y has been studied previously in the context of RS models in 


in the 5d picture and in 42 using a Kaluza-Klein mode decomposition. 42 maintains its 
focus on the decay B —)■ K* gT. Both works consider only the dominant effects of 5D 


41 


penguin diagrams and neglect the so-called wrong-chirality Higgs couplings terms 28 43 


This is equivalent to an RS model with a naively brane-localised Higgs that does not 
arise from a well-dehned limiting procedure. 

In this letter we want to consider the case of a bulk Higgs held. This scenario is quite 
general as it requires us to take into account both Higgs and KK Higgs contributions. In 
order to keep the advantages of the original setup, we still impose that the bulk Higgs 
is strongly IR localised. Following the construction of for the bulk Higgs gives a 5D 
wave function for the Higgs vacuum expectation value (vev) of the form 


viz) = 


+ /^) p/2y/3+ 

V 1 - e2+2/3 


K, ^/3+2 
USM ^ • 


(4) 


Here e = T/k, usm = n ~ 246 GeV and /3 is a parameter related to the 5D mass of the 
Higgs scalar. The typical width of the zero-mode prohle is determined by l/(/3T); the 
limit /? —)■ oo leads to a maximally localised ’bulk’ Higgs. For our subsequent analysis 
we always tacitly assume that this limit has been takenj^ We will focus on two types 
of RS models: the minimal model with the same gauge and fermion multiplets as the 
Standard Model and the custodially protected model with an extended matter 

and gauge sector (see 49 for details on the specihc setup). 

and we refer the 


This setup was also used in our work of lepton flavour violation 30 


reader to it and to 31 for explicit expressions for the 5D action and associated Feynman 


rules. For the study of the 6 —)■ sy transition we can directly transfer the results of 30 
to the quark sector. For simplicity, we only consider the effects of the strong interaction 
and the Higgs boson. Electroweak effects could be included in full analogy to the existing 
computation of flavour violation in the lepton sector, however, their inclusion will not 
lead to a fundamentally different phenomenology. Since we focus on QCD effects, we 
do not investigate the phenomenologically interesting decay B —)■ it receives 

tree-level contributions from four-fermion operators with both quark and lepton fields, 
which cannot be generated by gluon exchanges. 

The general strategy of the calculation then follows 


31 . We start with a fully 5D 


theory and integrate out the compact hfth dimension by matching onto an effective 


^See for details on how the limit has to be taken if 5D loops with a (dimensional) regulator are 
involved. 
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Lagrangian at the KK scale T. We will only consider operators of at most dimension 
six and the corresponding effective Lagrangian is the renowned Buchmiiller-Wyler La¬ 
grangian 50 . This step is presented in section 

We then transition from the symmetric phase to the broken electroweak phase. The 
Wilson coefficients of the resulting operators are subsequently evolved from the high 
scale T down to the typical scale of the process h —)■ sy, fib- This is discussed in Secj^ 
The phenomenological implications of the resulting corrections to the coefficients in the 
weak Hamiltonian are shown in section HI We conclude in SeclH 


2 Matching at the scale T 


A starting point for a completely general analysis of flavour-violating processes in BSM 
models is the Buchmiiller-Wyler Lagrangian [^. The new heavy degrees of freedom 
have been removed by matching onto the effective (dimension-six) Lagrangian. This will 
capture the dominant effects of any new physics model and only SM helds and dynamics 
are needed in any subsequent analysis. The price for taming a BSM model in this way 
is encoded in the (potentially) up to 2499 Wilson coefficient^ Each of these has to be 
determined by integrating out heavy degrees of freedom above the matching scale. 

Here we are only interested in the dominant contribution to 6 —)■ s transitions in a 
specihc class of RS models. That is, we only consider the flavour-changing transitions 
that are mediated by KK gluons and the (KK) Higgs. This greatly limits the number 
of operators that have to be considered. It is then convenient to consider the following 
effective Lagrangian at the KK scale fixK = T. 


pdimb — 


+ h.c. 


kDUi 




( 5 ) 


where we dropped operators that either will not contribute to leading logarithmic (LL) 
accuracy to 6 —)■ sy or are generated by exchange of SU(2), U(l) gauge bosons. Qi 
corresponds to a quark doublet of with generation index i; D and U are down- and up- 
type singlets. G and F are gluonic and electromagnetic held strength tensor, respectively; 

is a generator of SU{3) in the fundamental representation. Note that the Lagrangian 
is dehned in the unbroken electroweak phase and all quarks are still massless. Hence 
the indices i,j = 1,2,3 are not commensurate with e.g. up, charm or top. The ellipses 


all possible flavour structure are counted 


58 
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indicate a sizeable set of omitted operators that either cannot be generated via QCD 
effects or whose contribution to 6 —)■ S 7 is suppressed. 

In writing ([^ we tacitly assumed that we are in a flavour basis where the 5D fermion 
mass matrix is diagonal. Furthermore, ([^ already reflects the fact that we will need 
the coefficient of the electromagnetic dipole operator, i.e. instead of working with the 
field strength tensors of SU( 2 ) 2 , and U(l)y we only included the linear combination that 
will form the photon after EWSB. Using Fierz transformations it is possible to rewrite 
some of the operators in (|^ by removing the T“ ® T“ colour structure. This procedure 
is useful for a general analysis of flavour violation as one can use a minimal operator 
basis 


51 . For our simplified analysis this is not needed. 


The Wilson coefficients a and b will set the initial conditions for the RGE evolution 
from fixK to the electroweak scale fiEw ~ Mw ~ mtop where they will induce shifts in 
the coefficients of the well-known weak Hamiltonian. A subsequent evolution down to 
the scale Hb can then be performed in the standard way. 

Before moving on to the results for the matching calculation let us briefly review the 
parameters of the RS model that are relevant to our analysis. As for any BSM study of 
flavour the Yukawa matrices are of crucial importance. An RS Lagrangian incorporates 
two 5D dimensionless Yukawa matrices, Yu and Y^, corresponding to the couplings of the 
Higgs to up- and down-type SU(2)i singlets. We always impose that these matrices are 
anarchic, that is, the matrix elements are roughly of 0{1) and have random phases. Fur¬ 
thermore, as already mentioned above, each 5D Lagrangian (independent of the presence 
of a custodial protection mechanism) contains a 5D mass M^. for each 5D fermion field 
ijji. In practice, it is convenient to work with dimensionless parameters = M.^Jk. 
Hence, we have in total nine 5D mass parameters: cg., C[/., c^i- with z = 1, 2, 3. In order 
to obtain a phenomenologically viable low-energy theory that reproduces not only the 
SM quark masses but also the CKM matrix the mass parameters cannot be completely 
unrelated. E.g. c^., the mass parameters for the down-type singlets are usually not too 
far from —0.5. See 25 for details on the relation of the various parameters for anarchic 
RS models. 


2.1 Gluon-mediated four-fermion operators 

The simplest way to match the 5D theory in AdSs onto the effective Lagrangian is using 
5D Feynman rules 52 . This method is well established, see e.g. 29,31,35 for various 
applications. In particular, it avoids dealing directly with KK sums at the price of a 
more complicated integral structure in loops diagram. However, for operators that can 
be generated by tree-level interactions in the 5D theory no such complications occur and 
the matching calculation is straightforward. 

A further simplification for tree-level matching comes from the fact that there is 
only a very limited number of 5D vertex integrals that can occur. In particular, for 
intermediate gluons there is only one independent structure. The four-fermion Wilson 
coefficients differ only by symmetry factors and 5D mass parameters. We can then use 
the more general results of 30 for the matching onto four-lepton operators. One only 
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needs adjust the couplings and gauge-group factors accordingly. 
We find e.g. for the Wilson coefficient of the operator 

l/T 1/T 

l/k l/k 

where the zero-mode subtracted 5D gluon propagator i 
^ O.x.s) 


( 6 ) 


is given by 31 


+ 2x^ ln(a;T) 

T 

+ V ln( 2 /T) + 2?/^ In ^ I + 0{q^)^ + {x ^ y) (7) 


and the 5D wave functions are 

/£V) = s'XS^) = 


where e = T/k. 

Up to terms suppressed by the ratio e one then obtains 

uQD 


( 8 ) 


^ + &i(cqJ + hi{-CD,) + &2(cq,, Cd^) 


(9) 


with 


r, kis 

bn = -- 


2 1 


( 10 ) 


,-2 c _ d -1 


4 ln(l/e)’ 

gs^ {5-2c){l-2c) 6^-1 

4 (3-2c)2 

(1 - 2c<3)(12c£,)(3 - Cl-F cd) 1 e' 

2 (3-2cq)(3 + 2cl)(2-cl + cl) "" e 1 - e^cg-i i _ ^-20,-1 • 

The Wilson coefficients of all other operators are related to bf-^. They only differ 
by symmetry factors that take into account the exchange of identical quarks and the 
potentially different external wave functions and In particular, one hnds 



6g« = ^bf-^icD, ^ -cqJ 

( 11 ) 

“CliJ 

= b'^^{cQ^ -)■ -Cf/J . 

( 12 ) 


5 




















Figure 1: Topologies of 5D one-loop diagrams that contribute to the matching onto a'!^ 
at order ag. The external boson in the diagrams in hrst line can be a gluon or a photon. 
Internal bosons lines represent a 5D gluon propagator. 


2.2 Dipole operators 

The determination of the dipole coefficients a"' and a® is much more involved. Follow¬ 
ing the calculation of 31 the 0(as) contribution to a'^ requires the computation of the 


diagrams shown in the upper row of figure The contribution to involves the same 
diagrams (with the external photon replaced by a gluon) and the additional non-abelian 
diagrams shown in the lower row of hgure [T] Since the determination of the electro¬ 
magnetic dipole operators for leptons requires all topologies (see 31,^) both and 


can be obtained from known results by rescaling each individual diagram with a simple 
factors. This also implies directly that the 5D gauge invariance check for the leptonic 
calculation 31,32 can be carried over to the case of diagrams with (KK) gluons. 


Let us consider an example: The first diagram in the first row of figure [T] with both 
the internal and the external boson gluons. The contribution to a® can be obtained from 
the known result for same diagram topology with an external photon and an internal 
hypercharge boson B. Starting from this result we set all fermion hypercharges Yf to 
2, trade the U(l) couplings g' for Qs and replace the global factor iQfC from the photon 
vertex with —^iggT^. All other diagrams can be determined analogously. 

The way the computation of the dipole operator coefficients in 


31 is set up, we 


need to include contributions to the dipole structure from one-loop diagrams with an 
insertion of a four-quark operator, see Figure These extra terms ensure that the 
Wilson coefficient is scheme independent. This otherwise occurring scheme dependence 
is a well-known fact in flavour physics, see e.g. (60, 61 


By adding the contribution 
of the four-quark operators we can work with a scheme independent “effective dipole 
coefficient” analogous to the construction of (^. 


Due to the required chiral structure only four-quark operators that involve both 
doublet Q and singlets D can contribute: Qi'y^T'^UiDj'yDj. Up to a trivial colour 
factor this additional contribution is then again completely analogous to the one in the 
lepton case and we refer to for a detailed calculation. 
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Figure 2: Diagram topology that cancels the residual scheme dependence in . The 
box denotes an insertion of a four-fermion operator. 


2.2.1 Higgs contributions 

It is well known that loop diagrams with internal Higgs exchanges lead to a contribution 
to the dimension-six dipole operators that depends on products of three Yukawa matrices 
39 . This contribution can be sizeable and is an important source of flavour violation 


30 


28| . Therefore it is important to consider this Higgs contribution alongside the previously 
discussed gauge-contribution. 

For a bulk Higgs we further need to consider the effect of its KK excitations. The 
mass of the first few Higgs KK states is roughly proportional to the inverse width of 
the corresponding zero-mode 40,44,45 . Nonetheless these modes do not necessarily 


decouple even for a strongly localised zero-mode. This non-decoupling was first shown 


m 


40 ; the typical impact of the Higgs KK tower is comparable to the effect of the zero¬ 


mode alone and therefore non-negligible for the determination of the dipole operator 
coefficient. 

Let us first consider the effect of the zero-mode Higgs only. We can partially use 


the results of 30 for the leptonic dimension-six dipole operator to construct the corre¬ 
sponding result in the quark sector. Again we only need to replace U(l) charges and 
add SU(3) colour factors as appropriate. For diagrams where a Higgs is emitted from 
an external leg and not from the loop (see the diagram in £gure|^for an example), one 
further has to distinguish two different contributions: those where the external quark 
propagator propagates KK modes and so-called off-shell terms that arise if the external 
propagator is a mass-less zero-mode, but the 1/p^ pole in the propagator is cancelled 


by powers of p in the numerator, see 31 for a detailed discussion. The latter terms are 


basically irrelevant for leptons as they are effectively suppressed by a SM lepton Yukawa 
coupling. They may however play a role in the quark sector due to the large top Yukawa 
coupling and we include these terms in _ 

It is convenient to use the definition + iQfcA^ + igsT^G^ for the SM 

covariant derivative with being photon and gluon held; e is the charge of the 

positron. This definition then coincides with the choice usually employed in studies of 
the 6 —)■ sy transition, see e.g. 46,55 , and makes the expressions in the subsequent 


sections consistent with the standard literature. 
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In the minmal RS model we then hnd 


^ij\Higgs~ 19271'^ 2k^ Qu)Fq QdFd+{2Qe Qu)Fu) 

- (2Q^ + ««- Q‘) /a (i/r)inyjni«4’(vr) (w) 

“«I Hi,,, " I92ir2 F » ^ ^ 

“ ® /«(i/r)inn'ni« sihi/o. («) 

where the Fg are abbreviations for 

F, = /,g\i/T)ini«F(-c,j|Fi],i4"\i/T)=ini»jj‘“\i/r^ 

Fq = /r(i/r)|ni«9i“’(i/r)"[yi]t*F(caJ|n]ijg<">(i/r) 

f". = I^QAVT)[YiUF{-c,d[Y}A,,SZm^^^ 

Ft, = /Sh/r)|n|.tFT,(c,JlFj]^i4”’(l/r)2y4g<”'(l/T), (15) 

with 


k^ (3 - 2c) + (1 + 2c)e^'=-2 - (3 - 2c)(1 + 2c)e2‘=-i - (1 - 

(l + 2c)(3-2c)(l-e2'=-i)2 


Ft,{c) 7=i 


k4 I _ gl-2c 

1 - 2c 


(16) 


In writing the expression for Ft^{c) we assume that the mass parameter c is not too 
far from —0.5, which is realised for all parameter points that reproduce the low energy 
parameters of the SM. 

In the custodially protected model the Higgs contribution to the dipole is given by 


tiHiggs “ 19271^ 1:“^ 21:« Qu){Fq + Ft,) QdFd + i'^Q, Qu)F„) 

- W.I + 2Q. -2Qg) 4“'(i/r)inKin]ig9i“>(i/r) (i?) 

tlHigg. “ 192x2 ^ + 2^7-3 + Fd + F,) 

/ohvoinyni.g gfhim. (is) 


The terms in (13), ( [M] ), ( [l7| and ( p)^ with factors of Fg, q 
the off-shell contributions. 


T^,Q,u,d, correspond to 















Co.int 


Figure 3: Illustration the additional dependence of the 5D mass parameter dependence 
of the KK Higgs contribution for the diagram on the right. cq.„^, are the 5D masses 
of the internal doublet propagator and the external doublet zero-mode. For leptons 
generally only a small region in the upper right corner (cq ~ 0.4 —0.7) would be required. 


As already mentioned we also need to take into account the effect of Higgs KK modes. 
In we absorbed the effect of the KK bosons in global factors called i?*. These were 
assumed to be roughly independent of the 5D mass parameters and therefore allowed 
for compact analytic expressions. Nevertheless there is a nontrivial dependence of the 
KK contribution on the 5D mass parameters; in particular for diagrams with a Higgs 
emission from an external line. In the lepton sector this effect is quite small especially 
when compared to the sizeable numerical uncertainties; we therefore neglected it in 30 


In the quark sector the wide range of 5D masses leads to more noticeable effects; since 
we can only determine these numerically we do not give an explicit expression. To give 
an idea of the potential size: the left panel in £gure[^ shows the additional effect of the 
mass dependence (without numerical uncertainties) for the diagram shown on the right 
of the same hgure. One can see that the effect is indeed of the order a few percent for 
leptons, but can potentially be of 0{1) for quarks. It is therefore not feasible to use a 
simple analytic approximation as was done in the lepton sector. 

Furthermore, we need to include KK Higgs corrections to the off-shell contributions 
to the Wilson coefficients. Again these terms are not necessarily suppressed in the quark 
sector, as the third generation Yukawa couplings are sizeable. However, we can only 
determine this contribution analytically for the Higgs zero-mode and not for the Higgs 
KK modes; it is only accessible numerically, but is quite small, only about 25% of the 
corresponding zero-mode effect. 

We therefore treat the effect of whole Higgs KK modes similarly to how the gauge- 
contribution is handled. Here we only remark that the total effect of the KK modes is 


smaller than the effect of the Higgs zero-mode but not parametrically so, see also 30,40 
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2.3 Beyond QCD 


We only considered contributions to the Wilson coefficients that proportional to or, 
in the case of dipole operators, enhanced by 5D Yukawa couplings. Obviously, exchange 
of hypercharge bosons and SU(2) bosons will also generate four-fermion operators, con¬ 
tribute to both dipoles and give rise to operators of the schematic form 
The latter class of operators will contribute to e.g. flavour-changing Z couplings. 

The U(l) gauge coupling at a scale of 1 TeV is roughly a; 7 (i) ~ 0.01. The SU(2)i 
coupling is signihcantly larger with asu{ 2 ){f^ = 1 TeV) ~ 0.032, but still smaller than 
as{fi = 1 TeV) = 0.09. The fact that the weak coupling is only about a factor of three 
smaller than the strong coupling may warrant including weak effects in the matching 
calculation. Including the effect of the other gauge bosons is not a principle problem; 
their contribution to the four-fermion coefficients as well as the dipole coefficients can 
directly be obtained from results for leptonic dipoles in the literature, see (^ . 

A further effect that would have be taken into account when considering weak correc¬ 
tions is the modihcation of of SM parameters and relations that have been utilised in the 
SM computation. In particular the relation of and the W mass, that is frequently 
used when rewriting the SM expressions is affected by higher-dimensional operators 


see 


53 for the general case and 54 for the a discussion within the RS model). 


It should be noted that KK Higgses do not give rise to relevant contributions to the 
four-fermion operators if the Higgs zero-mode is strongly localised towards the IR brane, 
which we always assume. An exchange of a SM Higgs can give a contribution to the 
four-fermion operators. But only in a second matching step at the intermediate scale 
Hint ~ Mw when the Higgs degrees of freedom would be removed. In this case the flavour¬ 
changing Higgs coupling arise from dimension-six operators of the form Qi^Dj (see 
e.g. 1^). However, even then the contribution will be suppressed by an additional SM 
b-quark Yukawa coupling. We therefore ignore these contributions. 


3 Running to the low scale 


The typical energy release in a decay of the type B —)■ Y^y is of the order of the h quark 
mass and a typical scale choice is thus fib = Mb/2 ^ 2.6 GeV. From the Standard 


Model calculation of 5 —)■ sy in the framework of the weak effective Hamiltonian, see 55 
for an overview, it is known that the RGE evolution from the weak scale hw ~ 


down to Hb introduces sizeable operator mixing 56,57 


Our matching calculation was performed the scale hkk ~ T and QGD corrections 
are bound to be of importance. We then have two possible strategies: We can either 
evolve the terms in the dimension-six Lagrangian from the high scale hkk to the elec- 
troweak scale within the unbroken SM, then change to the broken phase and complete 
the evolution down to the scale Hb- The required anomalous dimensions for the hrst step 
can be found e.g. in 53,58,5^. Alternatively, we can work with the “broken” 


opera¬ 


tor basis already at the high scale and perform the evolution down to the low scale in 
one step (taking into account the top-mass threshold). The hrst approach is more in 
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the spirit of a matching onto a set of dimension-six operators. The second option has 
simpler “logistics” as we only need consider a single RGE. Both strategies are valid and 
ultimately must be equivalent in a situation where no additional dynamics between ^kk 
and fiw need to be taken into account. 

However, for the specihc process at hand the second option has the additional advan¬ 
tage that the structure of the required evolution equation has been studied in some detail 
in 63 . While 63 ultimately focusses on scenarios with e.g. a flavour-changing Z', their 


operator basis contains the full set of normal and colour-flipped four-quark operators. 
We therefore choose to follow this approach. 

Let us for clarity introduce the effective Hamiltonian at the high scale ^kk, that is 


used in 63 




AGf 

■71 


VtlVti, 


AC'i'^{flKK)Q7-y + ACsg{^KK)Q&g + ACj^{^kk)Q'7^ + AC'^g{^KK)Q'sg 


E E /\CI[A,B]{^fkk)Q\[A,B]+/\CI[A,B]{^fkk) Ql[A,B] 

A^B=L,R q=u,c,t,d,s,b 


ACt[A, B]{fiKK) Qi[A, B] + Ad^[A, H](/ixic) Qi[A, B] 


A,B=L,R 


where the operators are given by 


(19) 


Q77 = 



Qry = 



Q\[A,B] = (sc'y^PAbis) [qalgPEqa) 

Ql[A,B] = (Sa'y^PAba) {qalgPEqa) 


Qi[A,B] = (sc'y^PAdfi) {dfs'ygPRba) 

Q 2 [A,B] = {Sa'y^PAda) {dalgPEba) 

(20) 

with Pl/r = |(1 T 75 ) as usual and a, (3 
current-current and penguin operators 

are colour indices. Note that while the usual 

Ql = {Scl^PLCfi) [cn^gPiba) 

Q2 = {Sal^PLCa) gPlbn) 


Q3 = (SaYPLCa) i^ 3 lFPLqp) 

q=u,c,d,s,b 

Qa = {Sa^^PLCls) {qplgPLqo) 

q=u,c,d,s,b 


Q5 = (SaYPLCa) 

q=u,c,d,s,b 

Q% = [Sal^PlCfi) Y idplgPEqa) 

q=u,c,d,s,b 

(21) 

are not included in (19), they do enter the renormalisation group equations. This 

opera- 


tor basis is obviously non-minimal as e.g. Qi and Q^iL, L] are related via Fierz identities. 
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As we only consider the LO corrections due to new physics, this does not invalidate the 
RG analysis 55 . 

In total we have to consider 70 operators. Fortunately, there are only a few inde¬ 
pendent entries in the leading order (LO) anomalous dimension matrix. Most of which 
can be taken from 60,61 once the different operator normalisation has been taken into 
it|^ The remaining entries can be taken directly from 


account 


63 where the use of effec¬ 


tive, scheme-independent coefficients is implied. In the following we tacitly 

assume that C^/s refers to the effective quantity and forgo to display the superscript. 
We will not give the anomalous dimensions explicitly and refer to the original literature 
for details. 

With the anomalous dimensions at hand, the renormalisation group evolution equa¬ 
tion (RGE) 




( 22 ) 


can be solved in the standard way, provided the initial conditions at the high scale ^kk 
are known. As the anomalous dimension matrix 7 is sparse, a basis where the evolution 
is diagonal can be determined very efficiently. For the strong coupling constant we use 
as{Mz) = 0.1185 with decoupling of the top quark at rrit = 170 GeV. 

Once the evolution down to fib has been performed the result for the branching 


fraction of R —)■ Xs'y can be obtained using the formula 63,64 
Br {B —)■ 1 


Br (B ^ V7)||" + N 


+ + . (23) 


Here we use a minimum photon energy of R™”' = 1.6 GeV; the same as was used for the 
HFAG world average. Here the V is a non-perturbative correction 65 68 and we use 
N{E^ = 1.6 GeV) = 3.6 x IQ-^. 

Since we work in leading order in the new physics contribution, BSM effects only 
induce a shift in the Wilson coefficients 


d;’(w)^ c);'(w) +A7;((,i,), 

' L ' J SM ' 

The SM value of the dipole coefficients 

Cr-yifib) = -0.368 


k')/ 


k')/ 


(24) 


(25) 


can be taken from 15 . The primed coefficient is tiny as it is suppressed by nis/mb 
and can be neglected. 

For completeness we also give the formulae for the related process b —)• sg. It can be 
treated completely analogously; here the shifts in the coefficients Cgg(fib) are required. 


^In the corresponding operators Qis are only rescaled by a factor of 1/4 compared to their 
definition in (201,(21). The anomalous dimensions remain therefore the same. 
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The NLL SM prediction was determined in [^. The partial width for the process b ^ sg 
is given by 


T (5 ^ sg) 


a 




247r2 


|G 


fv:m 


|T>I 


I pbrems 
^ fin 


(26) 


The explicit expressions for 
is then obtained as 


and D can be found in 


69 


The branching fraction 


Br(6 —)■ sg) ~ 


r(6 —)■ sg) exp 

r(6^cez/e) ■ 


(27) 


With ^ 0.105 ± 0.005 being the experimental semi-leptonic branching fraction of 
the B-meson. In the SM one hnds 

Br^^(6 ^ sg) = (5.0 ± 1.0) x 10“^ . (28) 


The last missing piece for our analysis are the initial conditions for the RGE. That 
is, we need the Wilson coefficients AC in (19). 


Initial conditions 

The Wilson coefficient in at the high scale fixK can be obtained from the Wilson 

coefficients of the dimension-six operators in (|^. We need to rotate into the low-energy 
mass basis and replace the SM Higgs held (if present) by its vacuum expectation value 

v/V2: 


$ -)■ 



0 + \ 
{v + h + iC) J 


Qi —t Pl 


U^jUj\ 

uf/, J 


D^^ViPnd,. (29) 


We only take into account terms that contribute to the Wilson coefficients in (19) and 
drop all others. 

As an example, let us consider the ter m Di'f^T°‘DjUj in the dimension-six 

Lagrangian. Using the substitution rules (29) we hnd 


bV^ li°l,JYT-'Pnbqul,T-'Pnq^ = 

= S'fPftbq^q^PRq^ + SaYPRbf 

= 01" IR, R] + 0\'\R, R] (30) 


where a simple single sum over qu = u,c,t is implied. Here we dehned Id^bq^qu ~ 

general we will use the abbreviation 


oFF' _ \tdFZ fDF'lt hFP'FfP 
PaBCD — lAiV^ \CjOij ^iB^jD 


(31) 
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with the appropriate flavour rotation matrices R^^'\ 


Comparing (30) with (19), we obtain 


1 


72 




DU 


2'J'2 ^sbququ ’ 


72 


ACl-[R,R]{fiKK) = 


1 


7 


DU 


2N ^2^sbququ 

(32) 


The remaining four-quark operators can be related to operators in the weak Hamiltonian 
in the same fashion. For clarity, we have relayed the expressions for the Wilson coefficient 
of (19) to an Appendix. 

Similarly, we can obtain the effective dipole operator coefficients. Introducing the 
abbreviation we find 




V2 


KK) 


-- 

^GFV*Vtt 


IOtt^ 

em^T^ 

IOtt^ 


72 


+ E 


q=d,s,b 


Qq f nq Gf qf 

rj^2 >^qbsq 


mb 


72 


GsgifiRK) — 


erubT'^ 


IOtt^ 


a 


,711 
\sb 


a. 


AGFv:,Vtb^, 7, , 

-- ^Sg^f^KK) 


QsmbT^ 
IOtt^ 


72 


72 


gs mb T2 


a 


9 lt 

Jsfe 


V2 


E 

q=d^s^b 

- E 

q=d,Syb 

- E 


q=d,s,b 


Qq mq Cf oQD 

rp2 U sqqb 


mb 


m„ 




7, 


QD 


mn 


_ QD 


7, 


(33) 


All quantities on the right-hand side of (33) are implied to be evaluated at the scale 
^kk- The terms containing a /d-coefficient arise from the one-loop diagrams with an 
insertion of a four-fermion operator. They ensure that the (effective) coefficient ACj^ is 


scheme independent, see Sec. 2.2 


4 Phenomenology 

To see the potential effect of the additional contribution to G^'^ on the B —>■ Xgj decay 
we need to scan over the parameter space of the RS model. We will, as mentioned 
before, consider a minimal and a custodially protected RS model with an IR-localised 
bulk Higgs. The model parameters include the 5D masses of the fermions as well as the 
two Yukawa couplings and Y^. These parameters are not independent as we need 
to impose the condition that low-energy parameters of the SM are reproduced within 
uncertainties. We take into account the SM quark masses (at the scale T) and the CKM 
angles and phase; here we make use of the analytic approximations of [^. A further 
restriction is imposed by hand on the dimensionless Yukawa matrices as we require them 
to be anarchic. That is, the matrix elements all have roughly a common magnitude of 
0(1) and arbitrary phase. Similar to the analysis in we consider two samples of 
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Figure 4: Bt{B —)■ ^^ 7 ) as a function of the KK scale T. The blue (dark grey) points 
correspond to the data set with large Yukawas, Y^ax = 3. The orange (light grey) points 
correspond to Y^ax = 1/2. The horizontal lines indicate the experimental value of and 
uncertainty on the branching fraction. The left panel shows the result for the minimal 
RS model, the right panel for the custodially protected model. 


Yukawas: one with a maximum entry size of Ymax = 3 (representing the case of large 
Yukawa couplings) and one with an upper bound of Ymax = 1/2 (representing the case 
of small Yukawa couplings). 

The main result of our scan through the RS parameter space is shown in hgure It 
shows the branching fraction B —)■ Xg'y as a function of the KK scal^T for the minimal 
RS model (left panel) and the custodially protected model (right panel). The blue (dark 
grey) points correspond to Ymax = 3, the orange (light grey) points to Ymax = 1/2. The 
current experimental central value, see equation ([^, is represented by the solid horizontal 
line; the dashed lines indicate the uncertainty. 

We find that the branching fraction is, especially for small Yukawas, predominantly 
larger than in the SM. This is due to a sizeable contribution from that lacks an un¬ 
suppressed interference term with the SM contribution—its contribution to the branching 
fraction is always positive. In addition to that the contribution to is generally larger 
than the contribution to the unprimed dipole coefficient. The reason for this, as was 

is that the 5D profile of the doublet Q 3 (that very roughly 


41 


observed already in 

corresponds to the after EWSB) is typically larger than the prohles of the down- 
type singlets D near the IR brane; consequently the operator Q'^,^ oc {jR)aO'^''{bL)aF^y 
receives a larger BSM contribution. 

Only for the Ymax = 3 sample one can observe data points with a significantly reduced 
branching fraction compared to the SM. This is due to a destructive interference of 
and that can counteract the contribution due to C'r^^ if the Higgs contribution to 

is large. This effect is more pronounced in the custodially protected model 
where the additional fermion states enhance the dipole coefficient, cp. (13) and (17). 


For small Yukawas the phenomenology of minimal and custodially protected model is 
quite similar. This is to some extent a consequence of working only with QCD- and 


^Note that the mass of the first KK excitation of the gluon is roughly given by 2.5 x T 70 
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Figure 5: left panel: Effect of operator mixing on The histogram shows the 

distribution of without operator mixing relative to the full |AC' 7 ^(/if,)| with 

mixing (see text for details). The blue (dark grey) and orange (light grey) histogram 
corresponds to Ymax = 3 and Ymax = 1/2- right panel: Correlation of ACgg and 
in the custodially protected RS model for T = 4TeV. The triangle represents the SM 
values of and C^g] the dashed diagonal line indicates jACsgl = |AC 7 .^|. Same colour 
coding as in the left panel. 


Higgs-mediated contributions to the Wilson coefficient; QCD is treated the same in both 
models while the electroweak sector is extended and features additional bosonic modes. 
In the Ymax = 1/2 scenario the main distinction between the two models—the Higgs 
contribution— is suppressed. 

The smallness of the Higgs contribution for Ymax = 1/2 and the consequently smaller 
ACj^{fiKK) also make the inclusion of operator mixing mandatory. To see this we con¬ 
sider two quantities: the full /S.Cj^{pih) as obtained from the RGE (22) and /S.C'j^{pih)\naive 
which is also obtained via (22) but we set the Wilson coefficients of all four-fermion oper¬ 
ators at the high scale piKK to zero. We then consider the ratio /S.Ci-y{pif)\naivel 
The deviation of the ratio from one indicates the relative importance of the four-fermion 
operators for the 6 —)■ sy transition. Histograms of /S.C^-y{piif)\naivel are shown 
in the left panel of hgure For simplicity we only show the plot in the minimal model 
for T = 4TeV. For large Yukawas, Ymax = 3 in blue (dark grey), neglecting the con¬ 
tribution of from four-fermion operators leads on average to an increase of ^Cj^{pib) by 
5%. For a few Yukawa data sets the shift can be of the order of ±15%. In the case of 
small 5D Yukawa coupling (shown in orange) ignoring the four-fermion operator mixing 
basically always increases AC' 7 . This can lead to an overestimate of the BSM contribu¬ 
tion to the B —)■ As 7 branching fraction by up to 40%. Hence including the mixing is 
relevant and should not be neglected. This is of course quite general as FCNCs mediated 
by new, massive gauge bosons usually create simultaneous contributions to AC' 7 ..), and 
to the AC/ 2 [A, R] as is indicated by the need to include the four-fermion operators to 
obtain a scheme-independent result. 

For completeness we also show the correlation of AC 7 ^{p,b) and AC^^{p.b) in the 
right panel of hgure We see that on average the BSM contribution to is smaller 
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than the contribution to C^g as was also noted in This is more noticeable for the 
small Yukawa sample shown in orange (light grey). The two Wilson coefficients are 
then clearly correlated and one observes a ’’lower bound” on for a given value of 

ACg..),. However, with Y^ax = 3, it is straightforward to hnd parameter points where 
AC'y.y is much larger than the BSM contribution to Csg- The reason for this is the 
following: The zero-mode Higgs contribution to and are almost proportional to 
each other, see equations (|T3|)~(|g. However, the sizeable KK Higgs contribution has a 
more complicated structure; it contributes in a different way to and to a'^. This blurs 
the correlation. 

Finally, comparing with the experimental value for B —)■ Xg'y we hnd that for Ymax = 
1/2 the RS model parameter space is generally compatible with experimental data for 
T > 2 TeV. Since electroweak precision observables already put stricter bounds on 
the KK scale 25,26 , B —)■ does not give any new constraints on the KK scale. 

Nonetheless, sizeable corrections of about 5 — 10% are still possible. For large Yukawas 
the situation is much more intriguing, especially in the custodially protected model. As 
the large effects come almost exclusively from the Higgs exchange contribution to the 
dipole coefficients they are strongly dependent on the specihc form of the anarchic 
Yukawa matrices. It is difficult to deduce any hard bounds on the RS parameter space. 
However, the total BSM correction to the branching fraction can be quite substantial. 
Even for T ~ 5 TeV it is easy to hnd parameter points outside the current experimental 
limits. Consequently, the new Belle H searches would have the potential to discover 
the impact of KK states on R —)■ Xgj with masses well above 10 TeV. The search at 
the B-factory is therefore complementary to other powerful indirect search avenues like 
Higgs production and decay or dipole moments — experiments at vastly diherent energy 
scales. 


5 Conclusion 


We have studied the havour violating radiative transition 6 —)■ sy in RS models with an 
IR localised bulk Higgs. For simplicity, our analysis is restricted to QCD- and Higgs- 


mediated BSM ehects. We followed the strategy of 31 and matched the hve-dimensional 


RS model onto the SM ehective theory including dimension-six operators. Here we could 
make use of our recent results 30 for lepton-havour violation in the RS model. In par¬ 


ticular the complicated 5D loop integrals that determine the dimension-six quark dipole 
coefficients could be recovered from the electromagnetic dipole coefficient for leptons. 
This way we can include the effect of 5D loops with internal gauge, Higgs as well as KK 
Higgs bosons. 

After the transition to the broken electroweak phase, we used the results of 


63 to 


include the effect of operator mixing due to RGE evolution from the KK scale T to /r?, to 
LL accuracy. This is necessary as already in the SM the QCD corrections are sizeable and 
the dipole operator coefficient alone is not regularisation scheme independent. We hnd 
that for small Yukawa couplings, i.e., for small Higgs contributions to the dimension-six 
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dipoles, the mixing of additional four-fermion Wilson coefficients into can be sizeable 
and should not be neglected. We expect this to be true in any BSM model where dipole 
and four-fermion operators are generated via exchange of the same intermediate states. 

While our results for the Wilson coefficients are general, we assumed anarchic Yukawa 
couplings to study the phenomenology of the decay B —)• in both the minimal and 

the custodially protected RS model. We find that the additional contributions to the 
branching fraction can be sizeable for large Yukawas and moderate KK scales T. 

The strong sensitivity of the RS contribution to the specific form of the Yukawa 
matrices makes it challenging to directly constrain the parameter space of the model. 
Nonetheless, the decay is a useful tool that complements other powerful probes for the 
KK scale in the quark sector, like Higgs production/decay 1^,37,^. More importantly. 


for large 5D quark Yukawas there can be observable deviations of Br(R —)■ X^y) from 
its SM value even for masses of the hrst KK excitation of around 10 TeV. For small 5D 
Yukawas couplings {Ymax ~ 0.5) the impact of the RS model is mild; for KK scales that 
are not in conflict with electroweak precision measurements the B —)■ X^y branching 
fraction generally agrees with the current world average within uncertainties. In this 
case the aforementioned alternative search channels are more promising. 


Note added: While this work was in its final stage, 71 was published. 71 presents 


a detailed analysis of the 5 —?• sy transition in the minimal RS model with an exactly 
brane-localised Higgs. It is to our knowledge also the first computation of the RS contri¬ 
bution to dipole operators that does not rely on an expansion in the ratio of electroweak 
and KK scale. In addition to QCD and Higgs effects also electroweak effects are taken 
into account, but the model does, by construction, not involve Kaluza Klein Higgs con¬ 
tributions. 71 includes QCD operator mixing, but neglects the effect of the four-fermion 


operators. Since we consider the case of a localised bulk Higgs with KK modes, it is 
most useful to compare with the case of small Yukawa couplings; in this case the quite 
different Higgs sector does not play an all too dominant role. We then find RS corrections 
to Bt{B —)■ X^y) that are of similar but slightly smaller in size to those found 71 . This 


seems not unexpected as we neglect electroweak corrections to the dipole, but do in¬ 
clude mixing with dimension-six fermion operators, which tends to give rise to a slightly 
smaller ACj^ coefficient. 
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A Wilson coefficients of the extended electroweak 
Hamiltonian at the scale [Ikk 


In the following we collect the coefficients of the various four-fermion operators in (19). 
To this end we hrst map each operator in the dimension-six Lagrangian unto operators 
in the broken electroweak theory and extract the Wilson coefficients by comparing with 
(19). For brevity, let us hrst introduce the abbreviation V = 


bf QiYT^QiU,^;r-'Uj ^ = 

2 JY- (^sbququ Qu'y^P rQu T sbuu ^oTi^P l^ j3 (^Qu)P r(,(1u) a 

= Or|i. fl| + Ifti,.,. OML, R] (34) 

gives 

VAC?“[L.fl](,.,„0 = ^/3?41,. VACriL.R](,.KK)^-^l3Z,.- (35) 


= -YfS^OilR,R]+l}°°OflR,R] - ^ff^dilR,R]+l}SStOflR,R] (36) 
- jfffm OilR, R] + /3S£ R] - i/3r. OilR, R] + PS,Z 0\[R, R] (37) 


gives 


VACI[R,R]{^xkk) 

- —rdd 

rj-i2 ysbss 

VAC'2Ai?,/?](pxx) = - 

1 


V^Cl[R,R]{^iKK) 

- —RDD 

rj-12 sbbb 

V^CI[R,R]{ixkk) = ■ 

1 

WT2 

VACi[R,R]{iiKK) 

- —rdd 

rj-i 2 i^sbdd 

V^Ci[R,R]{^XKK) = - 

1 

iV,T2 

V/\Ci[R,R]{^lKK) 

- —rdd 

rjn2 ^sddb 

V^Ci[R,R]{^lKK) = - 

1 

iVcT2 


■A 


DD 

sbss 


rDD 

Psbbb 


nDD 

Psbdd 


nDD 

Dsddb 


( 38 ) 


b'^pQ.rT^Q.Qi'iRr^Q, 

= 0;|L, L] + OX[L, L] 
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( 39 ) 


gives 


Y&Oi[L,L] + flZ0i[L,L] - ^/3?iO|lL,L] +/3«« 0;[L,L] 
^fZ Oi\L, L] + fiZ, 01% L] - i/3f« 0\% L] + Z10\% L] 


oQQ 

• Z sbss 


VAC'i^[L,L](/i^^) 

_ ^ pQQ 

rj^2 Zsbss 

VAC'|[L,L](/i^^) = - 

1 

AeT2 

V/\C\[L,L]{iiKK) 

_ ^ pQQ 

rp2 Zsbbb 

VAC'2'[L,L](/iKx) = - 

1 

AcT2 

VAC'f[L,L](/i^^) 

_ ^ qQQ 

rj^2 Zsbdd 

V^Ci[L,L]{^^KK) = - 

1 

AeT2' 

VAaf[L,L](/i^^) 

_ ^ oQQ 
rj^2 Zsddb 

V^Ci[L,L\{^iKK) = - 

1 

A,T2' 


oQQ 

Zsbbb 


oQQ 

'Zsbdd 


VACML,L](ta) = ^);S«,. 


VACJ>[L,L](f,^K) = -^/3?0,. (40) 


Finally 

bZQiJ^T- 


+ \ZZO-\R,L] 

l^ZZ0i%R] + \l3Z0i%R] - i^fZ.,Oi\R,L] + \ZZ0i\R,L] 
l^fZ,di%R] + \ZM,di%R] - ^fZiOi\R%] + \ZSaOi\R,L] 
l^fZOi%R] + \ZZOl%R] - ^p%lOi[R%] + \Z.lOl\R%] 
^fZmO'i%R] + \l3Z0\%R] - ^/32SO$i>|fi,L| + \flZ,0'AR%] 

(41) 


gives 


V^Cr[R,L\{^^KK) = 


_rQD 

22^2 ZquQusb 

= ^ rQ^ 

2^2 ^ddsb 


Vl^Ct[R,L](p^% = 

Vl^Ci[R,L\{^XKK) 


_}_oQD 


97^2 ^sbdd 


■\-) A nd\ 
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VACt\R,L](p^KK) = :^_& 

vac; [L. «](,.,„,■) = ^/32f. 
VAC;|L,fl](M*-*-) = ^/?2S 
vac; [B.L] = ^/322 

VAC;|B,L1(m™-) = 5:^/?S 


VAC^[B, L](pkk) 
V^^|[L,R](|lKK) 

vac;[l,b](^ka-) 

VAC‘[L,fl](MAA) 

VAC|[fl.L](^AA) 

VAC2*lii,L](/JAA) 


1 

2iVcT2 

1 

2iVcT2 

1 

2iVeT2 

1 

2iVcT2 

1 

2iVcT2 

1 

2iVcT2 


/S 


QD 


/5 


Qr> 


f3 


QD 

sbss 




QD 

sbbb 


f3 


QD 

sssb 




QD 

bbsb 


(42) 
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